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c3 Abstract 

The Spin Hall Effect (SHE) is a promising way for transforming charge currents into spin currents 

^ in spintronic devices. Large values of the Spin Hall Angle, the characteristic parameter of the 

^ yield of this transformation, have been recently found in noble metals doped with nonmagnetic 

^ impurities. We show that this can be explained by resonant scattering off impurity states split by 

^ the spin-orbit interaction. We apply our calculation to the interpretation of experiments on copper 

o 

^ doped with 5d impurities and we describe the conditions to obtain the largest effects. 
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The Spin Hall Effect (SHE), first described by Dyakonov and Perel in 1971 [T], is a subject 
of intense research as it allows for the generation of spin currents in nonmagnetic conductors 
and the developments of spintronic devices built without ferromagnetic materials. The SHE 
is due to spin-orbit (S-0) interactions which deflect the spin up and spin down electrons 
of an electrical current in opposite directions. While the symmetry between spin up and 
spin down in a nonmagnetic material precludes charge accumulation on the edges of the 
conductor, i.e. no Hall voltage, there is a spin accumulation which can be exploited to 
generate a pure spin current. Alternatively a Hall voltage can be generated by injecting a 
spin-polarized current to break the symmetry, which is called the Inverse SHE effect. The 
SHE is associated with off-diagonal terms of the resistivity tensor having opposite signs for 
spin up and spin down electrons, respectively p^y and -pxyfor Sz = ±1/2. It can include an 
intrinsic contribution due to the effect of S-0 interactions on the wave functions of the pure 
material [21 [3] and an extrinsic one resulting from spin-orbit interactions on impurity or 
defect sites|H[3]. Two mechanisms can contribute to the extrinsic SHE, the skew scattering 
[1] and the scattering with side-jump [S]. 

When the SHE is used to produce a transverse spin current, the maximum yield of the 
transformation of a longitudinal charge current into a transverse spin current is related to 
the Spin Hall Angle (SHA), defined as = Pxy/Pxx where p^x is the diagonal term of 
the resistivity tensor, i.e., the conventional resistivity for spin a =^4, (±) channels. Conse- 
quently is the important parameter for practical applications in spintronics. Until 2007 
the largest values of obtained for pure materials, metals or semiconductors, had been 
obtained for Pt ^ 0.5% ) |6l [7]. The much larger value of 5% found in 2008 for Au [8] 
was surprising and has been ascribed to skew scattering by Fe or Pt impurities [9] . An even 
larger SHE 15%) was recently obtained by doping Au with Pt impurities [TO]. Actually 
this brings to mind the large values of of a few percent found thirty years ago [TTl [12] 
for the SHE induced by nonmagnetic 5d impurities in Cu, e.g. $// = 2.6%, for Cu doped 
with Ir. This large SHE, with a typical change of sign between the beginning (Lu) and the 
end (Ir) of the 5d series, was ascribed to resonant scattering on the impurity 5d states split 
by the S-0 interaction [H] . Recently measurements by Niimi et al [13] on Cu doped with Ir 
have confirmed the large value {^h ~ 1.5%) of the SHE induced by Ir in Cu and confirmed 
its skew scattering mechanism. Thus, impurity scattering appears as a promising way to 
obtain the most efficient transformation of charge currents into spin currents by SHE. This 
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has triggered the development of theoretical models of the SHE induced by impurities 

In this paper we present a calculation and discussion of the SHE induced by resonant 
scattering from impurity levels. For 5d impurities in Cu, we can explain the order of mag- 
nitude of the large SHE of the experiments and the change of sign between the beginning 
and the end of the 5d series. Whereas the recent papers of Gradhand et al. p!l] present 
calculations of only the skew scattering contribution, we calculate both the skew scattering 
and side jump terms. By comparing the spin Hall angles due to skew scattering and to 
side-jump, we can predict the threshold concentration at which the side-jump contribution 
becomes predominant and can generate very large effects. In contrast to the ab-initio cal- 
culations of Ref. [Hj, our calculation is performed in an analytical model which aims at a 
general description of the main features of the impurity-induced SHE and at a prediction of 
the best conditions for large effects. 

Our calculation is based on a partial wave analysis of the resonant scattering of free 
electrons from the j = 5/2 and j = 3/2 states of 5d impurities in a metal like Cu, as 
illustrated in the inset of Fig.l. From the splitting between the 5/2 and 3/2 levels, — 
E3/2 = 5Xd/2, where is the impurity 5d S-0 constant and by using the classical expression 
of the phase shift at energy E as a function of the resonant level energy Ej, ctn{rij) = 
{Ej—E)/A where A is the resonance width, we find to first order in Ac?/ A , At] = ''73/2— ^75/2 = 
5/2^sin^r72 where ■r]2 is the mean phase shift expressed as a function of the number of 
5d electrons on the impurity by Friedel's sum rule, ri2 = (3?73/2 + 2?75/2)/5 = ttZ^/ 10 jT3] . 
After expanding the states |j, rrij) in terms of |m, a) states and keeping only terms that will 
contribute to pxx and pxy ,we find the following expression for the scattering T-matrix (to 
first order in A^/A only non-spin-flip terms contribute). 
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where a = ±1 and n{eka) is the DOS for one direction of the spin. Note that interchanging 



k and k' in the first term in the bracket changes its sign; this is the signature of the anti- 
symmetric scattering. The second term is the usual symmetric term associated with charge 
scattering. 

From the antisymmetric part of the scattering probability Wantisymika — > k'a), associated 
with cross terms between the antisymmetric and symmetric parts of the T-matrix, we define 



iOskewikpa) by, 



^ Wantisym{^(y kV)5((k', Cx) =eU skew{kF(^)^ ■ k X Z, (2) 
k' 

where e is a unit vectors along the electric field E, z the spin quantization axis, tq the 
isotropic relaxation time and we have used the normal out-of-equilibrium distribution func- 
tion (7(k', cr) = —eTQVpG ■ k to arrive at this result. From the cross terms in the T— matrix 
between the / = 2 and / ± 1 terms and upon performing the integrals over spherical conduc- 
tion bands we find 

uJskew[kF(T) = -cr^^ — sin(2?72 - T]i) sm ?72 sm?7i, (3) 

where iVjis the number of impurities. 

The side jump contribution to the SHE enters when we write Hall current in the pres- 
ence of spin-orbit scattering [TB] as J = — e^lko-t'^k + i^a(k, c'"]/(k, a) where the term Ua 
is the anomalous velocity [arising from a side jump] attendant to electron fiow in systems 
with spin-orbit coupling [T7]. By using the distribution function /(k, cr) found from the 
linearized Boltzmann equation that accounts for antisymmetric (see Eq. 2) as well as sym- 
metric scattering and the appearance of side jumps [18], we find the current can be written 
as 



J(a =n) = (-^) roikFa)[vFk-Uaik,a)z x k)] 

X e ■ [vpk + {a;skew(^F, cr) + (^aikp, o-)}(z X k))]. 



(4) 



where we have written Wa(k, cr) = Ua{k,<j)k y< z . 

The transverse Hall current comes from terms proportional to z x e. When we consider 
spherical conduction bands and average over we find there are two transverse components 
in an = CTyx- The skew scattering one is 



CTskewi'y =ti) = ~3^^ 
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and the anomalous velocity or side jump contribution is. 
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FIG. 1: Skew scattering (squares) and side-jump (triangles) contributions to the Spin Hall Angle 
calculated from Eqs.9 and 13 as a function of the number of d electrons, Z^, for 5d impurities in 
Cu . The parameters are indicated in the text. The side-jump contribution is calculated for an 
impurity concentration of 2% . Inset: Density of States (DOS) of a 5d virtual bound state with 
S-0 splitting between j = 3/2 and j = 5/2 states. 

The normal conductivity for each spin channel is 



a 



-1 

N 



207ThNi 



2, sin^^a, (7) 

where kp is the momentum at the Femi level, and = \ntotai- By placing the expression 
for Uskewikpcr) , see Eq. 3, in the expression for (Tskew , Eq- 5, and dividing by aj^ we find to 
first order in (Ju/aj^ the skew scattering contribution to the Hall effect is. 



127riV,7iA. 



jSkewl^ =f I ^ ±) -_ _L 



sin(2?72 — rji) sin^ 772 sinr/i. 



The Hall angle from the skew scattering is 



A, 



^s^en^ = ±3/5^ sin(2r/2 - Vi) sinr/i. 



(8) 



(9) 



The contribution to the Hall effect from the side jump mechanism, i.e., the anomalous ve- 
locity, is found by repeating the calculation done for the Kondo-like rare-earth ions [16] , but 
this time using the T matrix for resonant d states, Eq.l. The expression for the anomalous 
velocity ^^(k, cr),Eq. (2.13) in Ref.l6 is. 
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Wa(k, a) = ^[Re VkTk^.kCT + Pi r ReT^^,^,^,Vk'Tkv/,k(T 

^ {£ka-£k'a') ' 

k/CT 



- TT ^ 5(efc<7 - efc'a')ImT,^^_,^,^, Vk'Tk'^'.ka]- (10) 

k/ff' 

Only the last term contributes to the Hall effect [16] and by using Eq. 1, we find 

Walk, cr) = cr — - — — — — — cos(3r72 - siu 772 sm T^ik X z. (11) 
iTn[eFjnkF A A 

By placing this expression in Eq. 6, and dividing by aj^ we find the anomalous velocity 
contribution to the Hall effect is, 

p~(a=ti^±) = ^ — — cos(3r/2-r7i)sm^r/2smr/i, (12) 

where c is the impurity concentration, and z = "'j^*"' = ^ , i.e., the number of conduction 
electrons per lattice site. Finally, the Hall angle from the side jump is, 

^anom ^ ^IQ /^^^^ ^^^^^^^ _ sj^S (^3) 

Similar calculations can be performed in the presence of crystal field. With completely 
crystal field split t2g and Cg states, for example, the prefactors of Eqs. (9) and (13) for the 
t2g states are multiplied by | and | respectively, and 772 is replaced by rjt^^ = ^Zt^^ , where 
Zt2gis the number of electrons in the t2g states. 

We begin the discussion of our results by a glance at the expressions of the Spin Hall 
Angle (SHA) for skew scattering and side-jump, respectively Eqs. (9) and (13), for the case 
without crystal field splittings. ^I^^'" is proportional to ^, and <|)^°''" ~ Large 
effects are thus expected for narrow resonances when the S-0 splitting induces significant 
differences in the scattering on 5/2 and 3/2 states. In the corresponding expressions for 
the intrinsic contribution to the SHA |2], |3] the denominator A is replaced by an energy of 
the order of the band width, therefore extrinsic effects due to resonant scattering should be 
generally larger in the usual case where the width of the resonance is smaller than the band 
width. 

The second important feature in the expressions for the Hall angle, arising from the 
symmetry rules for the SHE, is the interplay between the asymmetric scattering amplitude 
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in the channel / and the symmetric amphtudes in the channels / ± 1. It follows that the 
Spin Hall angle, Eqs. (9) and (13), depends not only on the phase shift 772 in the resonant 
channel (/ = 2) but also on the phase shift t]i in the non-resonant channel / = 1 (we 
have neglected the phase shift in the channel with / = 3). As the scattering in a non- 
resonant channel is generally weaker than in a resonant one, this selection of cross terms 
between different spherical harmonics (rarely described in theoretical papers) contributes to 
the general smallness of the SHA. 

We now focus on the skew scattering. If one supposes, as generally admitted, that 
the main contribution to the scattering by 5d impurities in noble metals comes from the 
resonance on their 5d states [15], rj2 is much larger then rji and, in first approximation, 
is proportional to sin 2772; see Eq. (9). As ri2 = sin 27^2 changes sign from positive to 
negative between the beginning and end of the 5d series as shown in Fig.l. This change 
arises from the difference in sign of the asymmetric resonant scattering on 5/2 and 3/2 
states. This agrees with the observed change of sign for the skew scattering SHE induced by 
5d impurities in Cu [III [12] , $h = -1-2% for Lu impurities {Zd[Lu] = 1) and = 2.6% 
for Ir impurities {Zd[Ir] = 8); the positive SHE for Ir in Cu has been confirmed by recent 
experiments [13j. A similar change in sign with Zd is observed also for pure 5d metals, 
which suggests a similar explanation based on the relative position of the 5d states, split 
by the S-0 interaction, with respect to the Fermi level. However,the crystal field splitting 
should also be taken into account for a precise prediction of the variation in the 5d series. 
As we have summarized after Eq. 13 the crystal field splitting between t2gSind Cg leads to a 
variation of the SHE as sin —j^ with a change of sign for Zt^g = 3. 

Now, we proceed to a quantitative discussion of the skew scattering Hall angle predicted 
by Eq. (9). First we discuss Cu doped with Ir for which different types of experiments have 
shown a predominant contribution from skew scattering with reasonably consistent values of 
= 2.6% in Ref.[n] and = 1-5% in Ref.|I3j . Typical values of A for 5d impurities 
in noble metals are close to 0.5eV from both experiments [I9] and ab-initio calculations 
P [20]. With A = O.beV, Xd ^ 0.25eV [21], and Zd = 8, the mean experimental value of 
the SHA for Cuir, = 2.05%, is obtained by introducing rji = —4.3° in Eq. (9). The 
decomposition of Eq. (9) into two factors, 3Ad^M2r?2-m) _ _o.277 and sinrji = —0.075, shows 
that the interference between the resonant and nonresonant channels induces a significant 
reduction. The calculation for Ir in Cu can be extended to other 5d impurities. With the 
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same values of A, 771, and using the S-0 constants for the 5d series [2T] and ■r]2 = one 
obtains the wavy variation of as a function of Zd shown in Fig.l. 

In contrast to the skew scattering contribution to the SHA, the side jump one is pro- 
portional to the impurity concentration c. The side-jump SHA for c = 2%, calculated with 
Ep = 7eV for Cu and the values of the parameters A^, A, rji already used for skew scat- 
tering, is compared in Fig.l with the skew scattering one. For impurities at the beginning 
and the end of the 5d series {Lu, Hf, Ir, Pt) the side jump contribution at c = 2% is much 
smaller than the skew scattering one. It is expected to remain smaller even at concentrations 
around 10%. This is in agreement with the results of a constant SHA up to c = 12% for 
Ir in Cu |T3]. On the other hand, for impurities in the middle of the series, like W,Ta or 
Os, concentrations as small as 2% yield side jump and skew scattering contributions of the 
same order of magnitude. For these impurities, very large contributions from the side-jump 
($^°™ > ]^g^ ^ expected for concentrations of the order of 10%. This can be compared 
to the situation of the Anomalous Hall Effect of Gd doped with Lu impurities [22] in which 
the side-jump contribution exceeds the skew scattering for concentrations above about 6% 
of Lu. 

The above discussion, for both the skew scattering and side-jump contributions, is altered 
when a crystal field splits the t2g and Cg states. According to the results summarized after 
Eq. (13), this introduces a change of sign not in the middle of the 5d series but at midway 
through the filling the t2g states at Zt^g = 3 , and a reduction by 3 and 5 in the amplitudes of 
the Hall angles. This can change the variation through the 5d series but not really the order 
of magnitude of the SHA's. For quantitative predictions only an ab-initio calculation of the 
scattering phase shifts can lead to realistic results. Our analytical calculation rather aims 
to predict the main features of what can be expected from 5d resonances and to identify the 
important parameters. 

Summarizing, large SHE effects induced by the resonant scattering from impurity states, 
here d levels, are expected from the combination of: i) a large S-0 coupling of the impurity 
states and a narrow resonance, which is the condition to obtain a large asymmetric scattering 
amplitude in the resonance channel I, and ii) a large symmetric scattering in the channels (Z± 
1). The second condition is not really fulfilled for 5d impurities in Cu so that skew scattering 
SHA's of only a few percent are expected in agreement with the existing experimental results. 
However, at least for some impurities {W,Ta,Os), large side jump effects can be expected 
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with exceeding 10% for c > 10%. We anticipate that similar features can also be found 
for impurities with p state resonances like Pb or Bi. Spin Hall angles above 10% would be 
extremely interesting candidates for generating spin currents without magnetic materials in 

spintronic devices. 

We thank Professor Elie Belorizky for helpful discussions on the projection of the scat- 
tering on crystal field states. 
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